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The  Line  Shape  for  CoUisionally  Broadened  Molecular 
Transitions:  A  Quantum  Theory  Satisfying 
the  Fluctuation  Dissipation  Theorem 


I.  INTRODUCTION 

The  dipole  autocorrelation  function  for  spectral  line  broadening  is  treated  in 

1  2 

a  quantum  theory  which  rigorously  satisfies  the  fluctuation  dissipation  theorem  ’ 
on  a  microscopic  level.  The  basic  approximation  in  the  theory  is  the  binary  colli¬ 
sion  approximation.  The  two-body  interaction  is  decomposed  into  a  part  which 
commutes  with  the  internal  coordinates  of  both  the  radiator  and  the  perturber,  and 
a  part  which  does  not.  The  theory,  as  developed,  is  appropriate  for  broadening 
mechanisms  for  which  the  non-commuting  terms  may  be  treated  within  the  frame¬ 
work  of  perturbation  theory,  while  the  commuting  term  is  to  be  treated  exactly. 

The  theory  gives,  at  long  times,  a  result  for  the  dipole  autocorrelation  function 
consistent  with  the  well -known  impact  approximation.  At  short  times,  an  auto¬ 
correlation  function  of  Gaussian  form,  with  a  renormalization  of  the  initial  state 
occupancy  is  obtained.  It  is  found  that  the  qualitative  features  discussed  above  are 
unaltered  in  higher  order  perturbation  theory.  The  results  are  consistent  with  the 
requirement  that  all  time  derivatives  of  the  autocorrelation  function  at  t  =  0  exist. 

This  further  satisfies  the  requirement  that  all  moments  of  the  lineshape  function 
in  the  frequency  domain  exist,  hence  that  the  lineshape  function  decays  "exponentially" 


(deceived  for  publication  24  September  1982) 
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In  order  to  pain  insight  into  the  influence  of  the  fluctuation-dissipation  theorem  on 
the  line  shape,  it  is  instructive  to  consider  the  effect  of  the  failure  to  satisfy  the 
FDT  on  the  far  wing  of  a  spectral  line.  Consider  the  impact  result  for  the  auto¬ 
correlation  function. 


0(t) 


Eh  ! i 


'yl'i 


if' 


fi 

e 


with  spectral  density 


♦  (a) )  H  i> . 

if  1 


>fi 


(u>  -  ulf.) 


(5a) 


(5b) 


where  p.  is  the  diagonal  element  of  the  density  matrix,  is  the  transition  moment 
describing  a  transition  from  the  i  to  f  state,  yfi  is  a  parameter  (the  halfwidth)  des¬ 
criptive  of  the  Markovian  collision  process,  and  uifi  (tiu>fi  =  (  f  -  t  d  is  the  transi¬ 
tion  frequency  (which  for  simplicity  includes  any  shift  term).  The  impact  result 
for  0<t)  clearly  does  not  satisfy  FDT  since  o(-t)  4  <J>(t  +  Mfi).  The  three  expressions 
for  the  imaginary  part  of  the  susceptibility  become 
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Thes.  three  impact  line  shapes  are  generally  referred  to  as  the  Van  Vleck- 
Weisskopf,  the  simple  I.orentz,  and  the  full  I.orcntz  respectively,  although  the 
front  dependence  on  fl  fiat  is  often  not  correctly  given.  If  we  consider  the  blue 
wings  for  which  fj-fiut  »1  and  u )  »  w)^.  we  obtain 


(2) 


(7a) 


x''(uO  =  L  Pi  lMifl 2  £ 

x>)  =  E  Pi  lpifl2  ~ 
\>>  =  E  Pi  lpifl2  { 


fi 
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This  clearly  demonstrates  that  an  autocorrelation  function  which  does  not  satisfy 
tlie  FUT,  for  example,  the  impact  autocorrelation  function,  gives  significantly 
different  results  for  the  far  line  wing  using  these  three  different  formulations. 
Additionally,  it  should  be  noted,  that  the  impact  autocorrelation  function  does  not 
satisfy  another  important  condition:  the  function  is  not  analytic  at  t  -  o  and  as  a 
consequence  the  higher  order  moments  in  the  frequency  domain  are  not  finite. 

An  example  of  the  type  of  problem  for  which  the  present  discussion  has  direct 

relevance  is  the  continuum  absorption  observed  in  the  self  broadened  water  vapor 
4  5 

spectrum.  ’  The  symmetrized  spectral  density  function,  <0 (a>)  4  p(-a>)a,  for 
this  case  is  shown  in  Figure  1.  The  water  vapor  continuum,  <oU>)  +  0(-*o)>con<i, 
utilized  in  this  computation  is  indicated  by  the  dotted  line  and  is  obtained  from  a 
line  by  line  calculation  from  which  the  contribution  within  25  cm  1  of  the  line  center 
has  been  excluded.  The  line  shape  is  taken  as  the  product  of  the  impact  line  pro¬ 
file,  Fq.  (5b),  with  an  exponential  type  function  containing  three  parameters  ad¬ 
justed  to  attain  agreement  with  the  continuum  observations.  The  continuum  based 
on  this  profile  provides  good  agreement  with  experiment  in  the  entire  spectral  re¬ 
gion  from  0-5000  cm  *,  including  the  three  atmospheric  "windows",  0-10  cm  ' , 
950-1250  cm  and  2000-2500  cm  '  as  shown  in  Figure  2.  For  comparison,  we 
show  the  continuum  resulting  from  the  impact  line  shape  alone  as  indicated  bv  the 
dashed  line  and  observe  that  <0(u,')  *  0( -u.')  >(.ont  is  loo  large  in  the  window  regions 
and  too  small  in  regions  of  strong  absorption,  it  is  clear  from  those  results  that  only 
a  proper  description  of  the  collisional  broadening  will  resolve  the  issue  of  the  role 
of  line  shape  vis-a-vis  other  mechanisms  in  contributing  to  the  continuum  absorp¬ 
tion  in  these  spectral  regions.  A  critical  test  of  the  theoretical  results  will  be  the 

extent  to  which  the  theory  describes  the  strong  temperature  dependence  of  the  ob- 

- 1  - 1  5 

served  continuum  absorption  in  the  1000  cm  and  2500  cm  regions.  ' 


4.  Clough,  S.  A.,  Kneizvs,  F.X.,  Davies,  K.  \V.  ,  Camacho,  l{.  ,  and  Tipping,  It.  II.  , 

(1980)  in  Atmospheric  Water  Vapor  edited  by  Deepak,  A.,  Uilkorson,  T.  D.  , 
Khunke,  I..  II.,  Academic,  I'ress,  New  York,  p.  25. 
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Figure  1.  Power  Spectral  Density.  <0(a>)  +  o(-u>)",  Self  Broadened  Water  Vapor. 
The  continuum  (...)  is  obtained  using  a  modified  impact  line  shape  (see  text) 

The  objective  of  the  present  theoretical  development  is  to  obtain  an  autocorrela¬ 
tion  function  which  satisfies  the  FDT,  treats  the  physics  of  collisions  and  is  com¬ 
putationally  tractable.  The  Hamiltonian  for  the  entire  system  may  be  separated 
into  two  parts. 


H  =  Ho  +  V  (8) 

in  which  Hq  describes  the  unperturbed  molecular  energies  including  intermolecular 
interactions  commutative  with  the  internal  coordinates  of  the  radiator  (including 
the  dipole  moment  operator  p)  and  the  perturber,  and  V  is  the  non-commutative 
part  of  the  interaction  potential.  The  matrix  elements  of  \  will  be  obtained  in  a 
representation  in  which  II  is  assumed  to  have  been  solved  exactly,  that  is,  in 
which  realistic  quantum  mechanical  trajectories  have  been  obtained.  The  theory 
is  developed  assuming  the  binary  collision  and  the  uncoupled  line  approximations. 
The  final  results  involve  a  second  order  perturbative  treatment  of  the  potential  \  . 


WAVENUMBER  (cm") 


Figure  2.  Power  Spectral  Density,  <  0(w,  +  p(-u>)> cont>  for  the  Self  Broadened 
Water  Vapor  Continuum.  Experimental  data  is  from  Burch,  Reference  5 


We  start  with  the  definition  for  the  dipole  moment  autocorrelation  function  for 
the  system,  oft ), 

0(t )  =  Tr  |M(o)  •  p(t)  p(H))  ,  (9) 

where  p(H)  is  the  canonical  density  operator,  and  JL(o)  is  the  dipole  moment 
operator.  For  /i(t)  given  by  the  Heisenberg  operator 

— ...  illt/fi  — ,  .  -iflt/fi 

M ft)  =  e  p(o)  e  (10) 

and  p(lt)  bv 

p(  H)  =  e'^“/Z  (11) 

with  Z  the  partition  sum,  we  obtain  for  o(t). 


/.v  ,-i.r  ,  ,  .  iiit/fi  — ,  . 

OU)  =  Z  I  r  [  /Ho)  r  4(0)  e 


Defining  a  time  development  operator,  L’(t)  through  the  relation 


-iHt /ft  -ill  t/li 


and  similarly  lor  complex  time,  t-i,j-h. 


-i  ll(t  -  i  A  h)/fi  -ill  t  ''-h  -  ,i  II 

'  it  '  f 


L  ( t - i  ,Mi> 


we  obtain 

Z  r  ill  t/fi  -ill  t/fi 

0(t)  =  Tr[iI(o>  •  U(t)  e  0  fj(o )e  0  p(HQ)  U(t-i/3f»)]  .  (15) 

At  this  stage  of  the  development,  we  sum  over  the  internal  radiator  states, 
designated  i  and  f,  suppressing  the  explicit  ni  dependence  for  simplicity  of  pre¬ 
sentation.  Uc  also  make  the  uncoupled  line  approximation,  although  this  approxi¬ 
mation  could  be  delayed  until  a  later  stage,  and  obtain  the  result 

Z  ,,  iuif.t 

0(t)  =  -77  £  'V  c:if(t)  e  (16) 

n 

whe  re 


C.f(t)  =  Tr[p(Ho)  U..  (t-i/lfi)  Uff(t)'|  . 


Here  H  designates  that  part  of  11^  not  dependent  on  internal  radiator  coordinates 
and  thus  the  trace  operation  is  to  be  taken  over  the  internal  perturber  states  and 
the  translational  states. 

The  binary  collision  approximation  assumes  that  the  correlation  function  for 
the  system,  C'^it),  is  expressible  as  a  product  of  correlation  coefficients,  c,j.(t), 
which  characterize  an  individual  radiator-perturber  pair. 


C.f(t)  =  |c.f(t>]  P 


where  N  is  the  number  of  perturbers  and 
P 

c.f(t)  =  Tr|ps<H()  U®  <t-i/3fi)  U®(t)t  ]  . 


(19) 


It  proves  convenient  to  define  the  operator  K^lt), 

Kf(t)  =  Tr[ps(lIo)  \V^f(t)T  j  (27) 

and  similarly  the  operators  K.  and  K.f  so  that  Cif(t)  is  expressible  as 

C.f(t)  =  exp  |  -  Np  [K.tt-id-fi)  +  Kf(t>  -  K.f(t, t-i/3fi>]  }  (28) 

and  pit)  as 

.  iu)  f  - 1  v 

0(t)  =  /'*  E  ()i  #i.f“  e  11  exp  |  -Np[K.(t-i/3fi)  4  Kf(t  >  -  K.f(t,  t-i/3-ft)]  J  (29) 


where  we  have  suppressed  the  superscript  s.  This  result  may  be  considered 
rigorous  within  the  binary  collision  and  uncoupled  line  approximations. 

The  evaluation  of  0 (t )  as  given  by  ICq.  (29)  remains  a  formidable  problem.  At 
this  stage  we  use  second  order  perturbation  theory  but  caution  that  the  procedure 
is  valid  only  if  the  matrix  elements  of  the  interaction  potential  V  are  small  in  the 
representation  in  which  I1q  is  diagonal.  We  proceed  by  evaluating  K^.(t )  through 
second  order  and  obtain 


Kf(t> 


OC  UG 

/  dkpk  /  d k-  (  ]/ 


ff 
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E 

j' 


(80) 


The  time  dependence  is  c  ontained  in  the  function  f  0(t)f 
1  ap 


V”  = 


~2  Id -cos  w^^t)  -  i  -  s.n  u)a(Jt|]  . 


(81) 


with 


a)  =  i,  -  t  i  .•(..  +  t,  -  t .  ,  .  (it 

aj3  f  f'  j  f  k  k1 

The  matrix  elements,  arc  obtained  in  the  representation  in  which 

IIQ  is  diagonal,  that  is,  the  representation  diagonal  in  the  internal  states  of  the 
radiator,  the  internal  states  of  the  perturber  and  the  translational  states  for  the 
isotropic  potential,  V  (  H).  For  the  dipole-dipole  interaction  appropriate  to  self- 
broadened  water  vapor  problem,  we  have  with  angular  integrations  and  m  sums 
suppressed. 
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where  p  and  p  are  unit  vectors  associated  with  the  dipole  moment  ol  «he  radiator 
and  perturber  respectively,  and  K  is  the  intermolecular  distance.  The  v.uvefunction. 

associated  with  translational  energy  <  is  obtained  from  the  solution  of  the 
Sehroedinger  equation  for  the  collisional  problem. 
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We  now  investigate  the  long  and  short  time  behavior  of  k^(t)  given  by  the  func¬ 
tion  f(t).  Making  use  of  the  identities 


lim 
t  -*  <* 


1  -  COS  ui  0t 
dtp 

- f - 


'  off 


1,1 


and 
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f  -»  0O 
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(25  a) 


(25b) 


with  Pr  denoting  the  principal  value,  we  can  obtain  the  long  time  behavior  of  Kf(t). 
For  the  real  part  we  have 


Re  Kf(t)  i  |t| 
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(36) 


and  for  the  imaginary  part. 


Im  Kj.(t)  °c  t 
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(37) 


This  result  for  long  times  is  consistent  with  the  time  dependence  obtained  from 

7 

other  formulations  of  the  impact  result  from  second  order  perturbation  theory: 
specifically  the  real  part  associated  with  the  width  varies  as  jtl  and  the  imaginary 
part  associated  with  the  shift  varies  as  t. 

7.  Anderson,  P.  W.  (1949)  Phvs.  Rev.  76:647. 
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l  or  sin. ill  times,  the  behavior  of  Kj.lt)  may  be  obtained  by  expanding  f(t),  and 
m  tin-  limit  as  t  -*  .1,  we  obtain 

y  X 

Kj.lt)  a  M  2t“/  elk  ,,  /  dk'  Z  X  2  .  (38) 

f  o  o  ff 

In  the  small  time  regime,  K.  (t )  is  a  real  quadratic-  I'unetion  of  time,  whereas 
K.(t-iph)  contains  a  contribution  from  the  FDT;  for  dipole-dipole  interactions,  Kjj. 
vanishes.2  The  real  part  of  the  correlation  function  in  this  limit  is,  nevertheless, 
>f  Gaussian  form. 

In  order  to  provide  additional  insight  into  the  functional  form  of  the  dipole 
moment  autocorrelation  function,  oft),  it  is  useful  to  consider  a  reference  function 
having  an  analytic  Fourier  transform.  We  choose  as  a  convenient  reference  the 
complex  function,  G(t),  for  which 

(  -at2  -aJtll 

K  .  =  ao  j  c  1  ©  e  J  (39a) 

and 

(  ^  -b2(th 

Im  G(t)  b}  |  te  ©  e  S  .  (39b) 

The  symbol  (x)  denotes  the  convolution  operation.  It  should  be  noted  that  a 

o 

is  chosen  such  that  G(o)  =  1  and  and  are  chosen  to  fit  the  short  and  long 

time  behavior  of  K  o (t )»  Eq.  (29).  Similarly,  for  Im  o(t),  b  and  b„  are  chosen 
o  o  2 

to  lit  the  short  and  long  time  behavior  of  Im  Q(t),  and  b;  is  assumed  equal  to 

Having  defined  the  reference  function,  we  now-  consider  the  difference,  A(t), 
between  the  dipole  autocorrelation  function,  0(t),  and  the  reference  function  G(t), 

A(t)  =  o(t)  -  G(t)  .  (40) 

This  procedure  has  two  purposes:  To  allow  consideration  of  the  results  of  the 
current  theory  on  a  scale  which  is  manageable,  and  to  facilitate  the  computation  of 
the  Fourier  transform  for  the  resultant  line  shape.  To  illustrate  the  functional 
properties  of  the  autocorrelation  function,  we  plot  in  Figure  3  a  typical  A(t)  calcu¬ 
lated  trom  an  approximate  0(t)  that  has  the  qualitatively  correct  time  dependence, 
l  or  a  more  complete  calculation  of  o(t)  we  expect  the  oscillations  to  be  less  pro¬ 
nounced  with  the  convergence  to  the  impact  limit  (KeA(t)=0)  to  occur  more  smoothly 
and  at  a  shorter  time.  The  rapid  variation  at  the  short  time  limit  is  expected  to 
remain,  reflecting  the  collision  dynamics  on  the  time  scale  commensurate  with  the 
duration  of  collision. 
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2»  c  t 


Figure  3.  Real  and  Imaginary  Parts  of  the 
Difference  Function,  A(t),  Calculated  From 
an  Approximate  0(t) 
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